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The Nambu- Jona-Lasinio model is supplemented by the quark interaction gen- 
erated by the one-gluon exchange. The employed gluon propagator exhibits 
the correct large momentum behavior of QCD, whereas the Landau-pole at 
low energies is screened. The emerging constituent quark model is one- loop 
renormalizable and interpolates between the phenomenologically successful 
Nambu- Jona-Lasinio model (modified by a transversal projector) at low en- 
ergies and perturbative QCD at high momenta. Consequently, the momen- 
tum dependence of the quark self-energy at high energy coincides with the 
prediction from perturbative QCD. The chiral phase transition is studied in 
dependence on the low energy four quark interaction strength in the Dyson- 
Schwinger equation approach. The critical exponents of the quark self-energy 
and the quark condensate are obtained. The latter exponent deviates from 
the NJL-result. Pion properties are addressed by means of the Bethe-Salpeter 
equation. The validity of the Gell-Mann-Oakes-Renner relation is verified. Fi- 
nally, we study the conditions under which the Nambu- Jona-Lasinio model is a 
decent approximation to our renormalizable theory as well as the shortcoming 
of the NJL-model due to its inherent non-renormalizability. 
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1 Introduction 



Most striking evidence for Quantumchromodynamics (QCD) as the right theory of 
strong interactions stems from high energy scattering experiments Jl], where exper- 
iments confirm the theoretical predictions. The remarkable property of asymptotic 
freedom implies that the running parameters of the theory, e.g. coupling constant 
and current masses, decrease with the inverse logarithm of the energy scale to the 
power of the corresponding anomalous dimension |]J , allowing a perturbative treat- 
ment. At low energies, the effective coupling constant becomes large and even di- 
verges at the so-called Landau pole which occurs at the typical energy scale of QCD 
Aqcd- It was recently argued ||, that the Landau pole is an artifact of the assump- 
tion of a perturbative vacuum, and that non-perturbative physics, which sets in at 
low energies scales, screens the Landau singularity. Since non-perturbative QCD 
is not yet under control, this implies that for a description of low-energy hadron 
physics one has to resort to effective models. 

Low-energy hadron physics is dominated by the approximate chiral symmetry 
embodied by QCD [BJ. In their pioneering work, Nambu and Jona-Lasinio pro- 
posed an effective theory with chirally invariant four fermion interaction [|J, and 
illuminated the mechanism of spontaneous breaking of chiral symmetry, which led 
to a successful description of the light pseudo-scalar mesons as Goldstone bosons. 
The simple handling of the model allows for a wide range of applications describing 



the light mesons |T0[ and baryons |JTl| . Recently, the model was also extended to 
describe heavy hadrons [|T2|]. 

Although the original version of Nambu and Jona-Lasinio was invented to de- 
scribe nucleons and their interactions, the re-interpretation of the model with quarks 
as fundamental particles gained further importance, when it was realized that the 
low-energy effective quark theory of QCD should be of NJL type. This idea is sup- 
ported in particular by the instanton picture || and the field strength approach to 
QCD @, 0]. In the context of strong coupling QED, it was argued that the so-called 
"irrelevant" four fermion interaction of the effective theory becomes relevant due to a 
large anomalous dimension || . It was further observed that a NJL- type fermion in- 
teraction naturally arises in the dual formulation of strongly coupled scalar QED ||. 
These field theoretical investigations indicate that an effective quark model should 
exhibit a NJL-type interaction in the low energy regime. 

Despite the wide spread success of the NJL-model, it suffers from two shortcom- 
ings, i.e. the absence of quark confinement, which e.g. allows for a decay of the heavy 
mesons into free quarks, and its non-renormalizability. The latter enforces that the 
interaction is cut off at some momentum scale, which defines the range of applicabil- 
ity of the model. As a consequence, processes involving large momentum structure, 
e.g. form-factors at high momenta and anomalous decays like 7r — > 77, cannot be 
adequately described by the model. Physics which involve the interaction at high 
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momentum can be properly described in the one-gluon exchange models [13, 114, [15 
which match the correct high energy behavior of QCD. 

The confining mechanism of QCD is still an unresolved problem. In certain 
super-symmetric Yang-Mills theories, it became apparent that the condensation of 
monopoles plays a key role jTjJ. Unfortunately, a concise understanding of confine- 
ment in QCD is not available at the moment. At the level of an effective quark 
theory, an extension of the NJL-model was proposed which realizes confinement by 
random background fields jl7|- In the context of the one-gluon exchange models, 
confinement is described by a quark propagator which vanishes due to infra-red 



singularities ||14|| , or which does not have poles corresponding to asymptotic quark 
states ]nj. In this paper, we will not address the question of confinement at the 
level of effective theories, but focus on the consequences of renormalizabilityQ. 

Here we propose an effective quark model which in the low energy regime co- 
incides with the NJL-model modified by a transversal projector (see equation @), 
and which matches the well known perturbative quark interaction at high energies 
at an energy scale Aj. We naturally choose Aj > Aqqd in order that the NJL-part 
of the interaction screens the unphysical Landau pole. Since the interaction strength 
of the NJL-interaction is fixed by the constraint to match the interaction induced 
by the one-gluon exchange at high energies, the only parameter of our model is Aj. 



Our model should not be mixed up with the gauged NJL-model |8|, [TJ, [T9| . In 
the latter, the NJL-interaction exists over the whole momentum range in addition to 
the gauge interactions. This model is non-perturbatively renormalizable and repre- 
sents an interesting alternative to the standard mechanism of electro- weak symmetry 
breaking [ED| . It is, however, not suitable in the context of the standard formula- 
tion of QCD. This is because its effective quark interaction does not reduce to the 
one-gluon exchange at high energies. 

The organization of the paper is as follows: In the next section, our model is 
defined, and its renormalization is discussed at the level of the Dyson-Schwinger 
equation. To this aim, the high-energy behavior of the quark self-energy is investi- 
gated by solving the Dyson-Schwinger equation numerically as well as analytically, 
using the Landau approximation in the latter case. Furthermore, the dependence 
of the chiral phase transition on the low energy four quark interaction strength is 
studied in some detail. In section [3], pion physics is addressed by means of the 
Bethe-Salpeter equation. The Gell-Mann-Oakes-Renner relation is established. Nu- 
merical results concerning pion properties are presented. Finally, we compare our 
model with the standard NJL-model in section |j. In particular, the shortcomings 
of the latter model due to its non-renormalizability are discussed. A summary and 
conclusions are left to the last section. Some explicit calculations are given in three 



1 At the moment, there is little empirical evidence that confinement is important in hadron- 
physics at energy scales below the would-be quark thresholds. 
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appendices. 



2 Model description and renormalization 



2.1 Model building 

For the description of hadrons we propose an effective one-gluon exchange model 
which is defined by the following generating functional (in Euclidean Space) 



s 



VqVq exp{- 



S 



d x q(x)(f)(x)q(x) 



d A x q(i$ + im)q + | / d*xd 4 y q{x)t a ^, Jb q{x)V'^ v {x - y)q(y)t b -i v q(y) 



where m is the bare current quark mass, t a , a = 1 ... 8 are the generators of SU(3) 
color group and 7^ are the Hermitian Dirac matrices. Furthermore, <p(x) is an 
external meson source. For m — 0, the model is chiral invariant. The interaction 
is diagonal in color space, V ab v = 5 ab D^ u , and consists of two parts, which, in 
momentum space, are separated by an energy scale A/, i.e. 

4ira(k 2 ) 



G0(A/- \k\) + 



k 2 



A, 
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k^ky 

~k 2 ~ 



(2) 



where 9{x) is the step function and | • | denotes the module at a Euclidean 4- vector. 
Choosing 

a{k2) = (33-2N f )\nk 2 /Al CD ' (3) 

where Nf is the number of quark flavors, the second term in (||]) represents for k > Aj 
the perturbative gluon propagator of QCD multiplied by the effective (running) 
coupling constant a(k 2 ). The coupling strength G of the NJL-term is chosen to 
match the gluonic interaction at the scale A/, i.e. 

4vra(A2) 



G 



A? 



(4) 



In order to study the spontaneous breakdown of chiral symmetry, we investigate 
the quark self-energy S(p). It can be obtained from the Dyson-Schwinger equation, 
which in ladder approximation reads 



Here 



ij> + m + - J ^34 lnS(k + p)-i v D^(k) = S \p) . 

Z(p 2 )^> + zE(p) 



(5) 
(6) 
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is the quark propagator, which is diagonal in color space, and we have used 
Sa,z ^ifith = I ^ik- K will turn out (see section |27|) that the loop integration in (^) is 
in fact ultra-violet finite, since the gluon propagator decreases like 1/k 2 ln(fc 2 / Aq CD ) 
for large values of k 2 , and the self-consistent solution £(fc 2 ) is asymptotically pro- 
portional to l/[ln k 2 //j, 2 ] dm , where d m = 4/9 is the anomalous mass dimension for 
three quark flavors. In addition, no finite renormalization constants are needed in 
(|j), because the freedom in the functions Z{p 2 ) and £Qo 2 ) is sufficient to accomplish 
for the standard renormalization conditions of the quark propagator. 

Using the ansatz (|6|), the Dyson- Schwinger equation can be cast into two equa- 
tions which determine the self-energy SQo 2 ) and the quark wave-function Z(p 2 ), 
i.e. 



Y( 2 \ AC I 24 (JO + K) ) . . 

U ' n m+ J\k\<Aj J^Y Z 2 ((p + k) 2 ){p + k) 2 + £ 2 ((p + A;) 2 ) M 
d 4 k £((p+fc) 2 ) Ana(k 2 ^ 



+ 4 



fc |>A, (2tt) 4 Z 2 ((p + k) 2 )(p + k) 2 + E 2 ((p + k) 2 



3 7| fc |< A/ (2vr) 4 Z 2 ((p + k) 2 )(p + k) 2 + S 2 



4 r _d*k Z((p + k) 2 ) ( 

+ 3 i|fc|>A / (2^)4 Z 2 {{p + k) 2 ){p + k) 2 + T 2 { 

2\ 




4na(k 

x ^ 

Note, that the k- integration never crosses the Landau singularity at \k\ = Aqcd due 
to our assumption A/ > Aqcd- 



2.2 Numerical results 

To pave the ground for later investigations, we first present a numerical solution 
of the Dyson-Schwinger equation in the chiral limit m = and discuss its relevant 
properties. The numerical method to solve the system of coupled differential equa- 
tions as well as the procedure to estimate the numerical errors are described in some 
detail in appendix 

The actual choice of the fundamental energy scale Aq C d depends on the reg- 
ularization scheme, which is used to extract the divergent parts. In the MS- 
scheme, which employs dimensional regularization, its typical value is Aqcd ~ 
100 . . . 200 MeV ffl, whereas in the one-gluon exchange models, which generically 
use a sharp 0(4)-invariant cutoff (and therefore a momentum-subtraction scheme), 
values up to 500 MeV are consistent with phenomenology |T3], |TJJ . In the case of 



the present model, it will turn out that the choice Aqcd — 500 MeV reproduces 
the experimental data for the static pion properties (see section |3|). For illustrative 
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Figure 1: The self-energy £(p 2 ) times 1 + p 2 / Aq CD in units of Aq C d and the wave 
function Z(p 2 ) as function of the momentum transfer p 2 (full lines). Also shown are 
the contributions from the NJL part of the interaction (long dashed) and from the 
gluon exchange term (dot-dashed). 



purposes, let us choose G = 1.11 16tt 2 / Aq Cd . It will turn out, that this choice for 
G is to small to provide a constituent quark mass of M ~ 310 MeV. 

The numerical result for the self-energy £(p 2 ) and the wave function Z(p 2 ) for 
three flavors (Nf = 3) is shown in figure |]. Although both integrals in (0) are equally 
important in order to calculate the self-consistent solution £(p 2 ), it is instructive to 
split S(p 2 ) into the contributions stemming from the NJL type interaction (first line 
in (|7|) ) and from the one-gluon exchange term after the self-consistent solution was 
obtained. These contributions are also plotted in figure [I]. As expected, the NJL 
term dominates at small momentum transfer, whereas the high momentum behavior 
is exclusively determined by the one-gluon exchange part of the interaction. 

The high-momentum behavior of the quark constituent mass function, i.e. 

MV) := gg , (9, 

is of particular interest, since it can be compared with the result from perturba- 
tive QCD. Including the leading operator-product corrections, the mass function is 
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Figure 2: The asymptotic behavior of the quark mass function M(p 2 ) (points) 
compared with the result of perturbative QCD (line). 



asymptotically (for large values of p 2 
2 , rn R (n) \\n(fi 2 /A 2 QCD ) ld: 



M{p z ) 



given 



An 2 d r . 



by m, 

, (qq)(/*) [ln(/i 2 /A 2 



QCD) 



-dm 



HVV^QCD, 



3p 2 



Hp 2 /a 2 qcd ) 



l-dr, 



+ 



(10) 



where ttir and (qq) is the renormalized current mass and the quark condensate, 
respectively, at a renormalization point /i, and d m is the anomalous dimension of 
the current mass. In perturbative QCD one finds 

1 2 4 

N f = 3. (11) 



33 - 2N 



f 



9 



for 



In the chiral limit, the mass function M(p 2 ) asymptotically decreases like 1/p 2 with 
logarithmic corrections. The same is true for the self-energy S(p 2 ), since the wave 
function Z(p 2 ) logarithmically approaches 1. Figure [| exhibits the asymptotic be- 
havior of the numerical result for M(p 2 ). The numerical data show a logarithmic 
decrease of p 2 M(p 2 ). The line corresponds to an anomalous dimension of the quark 
condensate of 1 — d m = 5/9, which agrees with the result of perturbative QCD JE]. 



7 



2.3 Analytical results 



The asymptotic behavior of the quark mass function ( |TUD is well known from per- 
turbative QCD studies augmented with OPE corrections ||21|| . In the following we 
will show that our model Dyson- Schwinger equation (^,|]) yields a quark self-energy 
XQo 2 ) which has the correct asymptotic (large p 2 ) behavior predicted by QCD. The 
asymptotic form of Dyson- Schwinger equations was studied in detail for many mod- 
els. For a review we refer to the book of Miransky |22 . 



In the case of our model, we firstly show that the NJL part of the integral equa- 
tion (0) does not affect the asymptotic behavior of £(p 2 ), as it was already observed 
in the numerical calculation (see figure |l|). In fact, since the momentum integration 
in the NJL-term is over a finite volume in momentum space, the integrand can be 
expanded in powers of k 2 /A 2 . Obviously, this contribution to S(p 2 ) is suppressed 
by powers of A 2 /p 2 , i.e. 

G Aj 



5V 



£(p 2 ) + . . . 



(12) 



8tt 2 Z 2 (p 2 )p 2 

implying that asymptotically (p 2 — ► oo) only the gluon exchange term contributes. 

In appendix |B| we show that in the Landau approximation^ and for non- vanishing 
renormalized current mass the coupled set of Dyson-Schwinger equations (^,||) has 
for p 2 — > oo the asymptotic solution 



K 



In 



A 2 



where 
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a 



d„ 



(13) 



(14) 



33 - 2N f 

and k is proportional to the renormalized current mass tur. Furthermore, the bare 
current mass runs with the cutoff as 



m(A uv ) = — , k = m R (fi) [ln(/i 2 / ' A 2 QCD 

[ m ly UV I ly QC 'D\ 

In the chiral limit (m# = 0), the asymptotic solution is given by 



P z 



In- 



'QCD 



(15) 



(16) 



where c is measure for the quark condensate (compare with equation (|I~0D). 

The main observation in this subsection is that the large momentum behavior 
of the solutions S(p 2 ), Z{p 2 ) of the coupled set (0,|) of differential equations agrees 
with the QCD predictions. 

2 This approximation is only applied in this subsection. The numerical studies in this paper 
treat the full Dyson-Schwinger equations 
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Figure 3: The quark mass function M(p 2 = 0) (crosses) and Z(p 2 = 0) (crosses) at 
zero momentum transfer as function of the coupling strength G of the low energy 
effective four quark interaction. Also shown is the fit (solid line) according fll8|). 



2.4 Chiral phase transition 

The results found in the previous subsection, i.e. a non-vanishing quark self-energy 
with the correct decrease at asymptotic large values of the momentum transfer, 
establishes that chiral symmetry is spontaneously broken for an appropriate choice 
of parameters (G > G c ). Here, we will report details of the chiral phase transition 
at the critical value G c of the effective four quark interaction. We have numerically 
investigated the quark self-energy S(p 2 = 0) and the wave function Z(p 2 = 0) at 
zero momentum transfer by solving the coupled set of integral equations (0) and 
(|]) for various values of G, the only parameter^ which is not fixed by perturbative 
QCD results. The results are depicted in figure [| For small values of G < G c , the 
chiral symmetry is not spontaneously broken and E(p 2 ) = is the only solution to 
equations (0JD- However, if the coupling exceeds the critical valuef] 

G c « (0.26 ± °- 011 ) x T?- > ( 17 ) 

^ ' A-QCD 

a non-trivial self-energy signals the spontaneous breakdown of chiral symmetry. For 
values G slightly above the critical coupling G c , our numerical results are well fitted 

3 Note that Aj is fixed by the constraint (f|). 

4 The procedure to obtain an estimate of the numerical error is discussed in appendix H. 
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Figure 4: The quark condensate x = ~ (qq) a t a renormalization point \x = 1 GeV 
as function of the low energy interaction strength G. The solid line in the right hand 
picture is the fit function fl20|). 



by 



M(p 2 = 0) 



0.43 A Q cd [pr 



5, 



0.47 



If e = G/G c — 1 measures the deviation of the coupling strength from its critical 
value, the behavior of the order parameter close to the phase transition is in general 
M ~ e s , where the critical exponent is characteristic of the underlying field theory 
In our case fllTf), a critical exponent of 1/2 is compatible with the numerical 



only p5 
data. 

Figure [3] also shows the wave function Z(p 2 = 0) as function of G. For G < G c , 
the wave function takes its canonical value Z(p 2 = 0) = 1. For values G larger than 
the critical coupling G c , the function Z(p 2 = 0) linearly increases with increasing 
strength G. It saturates, and becomes a constant independent of G for sufficiently 
large G. 

Finally, we present the quark condensate as a function of the interaction strength 
G. The quark condensate at a renormalization point fi is defined via the trace of 
the quark propagator S(k) (0) (see subsection |3.2|) , i.e. 



(w)0) 



In 



1 d n 



A 2 

1V QCD 



A 2 
In Auv 



A 2 

1V QCD 



fc2<Ar 



d A k 



trS{k). (19) 
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It was shown in |26| that this definition of the quark condensate is compatible with 
the definition via the asymptotic behavior of the mass function ([K]), which is used 
in the operator product expansion (see also subsection pT2|) . The numerical result 
of the quark condensate at a renormalization point of fi = 1 GeV is shown in figure 
|j. For values G close to G c , the numerical data is well represented by the function 

(qq)(ji = 1 GeV) « -0.041 A 3 QCD (JL - l)^ , with 5 X « 0.27 . (20) 

In the standard NJL-model, the quark condensate is directly proportional to the 
self-energy implying that this model yields a critical exponent of 1/2. This simple 
connection between the self-energy and the condensate is artificial, since the conden- 
sate is in general a subject of renormalization. In particular, the critical exponent 
of our model significantly deviates from 1/2. 



3 Static pion properties 
3.1 The Bethe-Salpeter equation 

In the constituent quark picture, the pion is described as a pseudo-scalar bound 
state of a quark anti-quark pair. This bound state is characterized by its vertex 
function P(p,q), which satisfies a Bethe-Salpeter (BS) equation. The momentum 
p thereby represents the four momentum of the pion, and q reflects the internal 
structure of the pion. In our case, the BS-equation is 

P(p, q) = I i 7 „ ta D A k - Q) S(k + p/2) P(p, k) S(k - p/2) t a i lv , (21) 

Jk 

where D^k — p) is the effective interaction (0) of our model. 

In the following, we will solve this equation ([H|) within a derivative expansion. 
This expansion was earlier obtained in references |£J and [23|. For later use, we 
briefly re-derive the expansion in our notation. 



For this purpose, first note that the BS-equation fl2"T| ) can be formally written as 
an eigenvalue equation, i.e. 

P p = [ K(E,p 2 ,k,q) P p (q) or: P p = K(X,p)P p , (22) 

Jg 

where the BS-vertex P p is represented as a vector, and the kernel K(E,p) depends 
on the self-energy S(A; 2 ) and the pion momentum p. For later convenience, we 
introduce a "scalar product" of two vectors a and b via a quark loop integral, i.e. 

(o|6):=tr J-0^S o (k-^a\k-^k + ^S o (k + ^b(k + ^k-^). (23) 
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As a warm-up exercise, we show that the pion BS-vertex is given in the chiral limit 
by 

Po(k 2 ) = ^ 75 Z Q (k 2 ) , (24) 

JTT 

where S is the self-energy in the chiral limit, and f n is the pion decay con- 
stant. The relation (p4|) is a particular case of the general relation between the 
BS-equation at zero energy and the Schwinger- Dyson equation. From (^|), one has 
that S {k)j 5 S (k) = ^(Zl(k 2 )k 2 + Eg (A; 2 )) -1 . A direct calculation then yields 

K(S ,p 2 = 0)P = Po, (25) 



where the Dyson-Schwinger equation (Rl) was used. Po (24) indeed describes a 
pseudo-scalar quark anti-quark bound state with zero mass, i.e m 2 = — p 2 = 0. 
The normalization of the BS-vertex is not fixed by the linear equation (P2"Q. It can 
be related to the pion-decay constant by assuming that the axial- vector quark anti- 
quark coupling is entirely dominated by the axial-vector pion coupling in the chiral 
limit |15| . 

We now consider the case of a small explicit breaking of the chiral symmetry 
by the current mass m^. In order to calculate the pion mass treating the explicit 
chiral symmetry breaking as a perturbation, we expand the BS-equation (|22|) to first 
orderQ in p 2 = —m 2 and ttlr, i.e. 

K(E ,p = 0) + K'(J: ,p = 0)p 2 + 5K](Po + 8P) = (Po + 5P), (26) 

where K' denotes the derivative of K with respect to p 2 , and 

SK := K(H,p = 0) — fC(S ,p = 0). The change in the pion mass in first order per- 
turbation theory obtained by sandwiching equation (^) with Pj. The contributions 
containing SP drop out as a consequence of the BS-equation in the chiral limit, and 
we are left with 

= -(P | #'(£„, 0)P ) ml + (P | 5KP ) . (27) 
In order to evaluate the second term at the right hand side of (p7|), one uses 
(K(S,0) - K(S ,0))P = 

^ ' (H ' i^t a D flu (k-q)[S(k)S(-k) - S (k)S (-k)}Z (k) l5 t a i lu 



U J (2vr) 4 
to obtain finally 



(p \skp) - 8Nc f d ' k s2f ^JKriMi^_ m 

(jPo|<5KPo) " IfJWji o( \z 2 k 2 + Y?)(Z 2 k 2 + Eg) ' (2 * 



5 Note that current algebra shows that is of order m^, which is also assumed in chiral 
perturbation theory. 
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The first term at the right hand side of fl27|) is more involved. This term explicitly 
reads 



(P o \K'(E o ,0)P Q ) = (Pi 



dK(E ,p 2 



dp 2 



p 2 =0 



P o) = p >< 

d 



(29) 



tr / t a i lv S(q) l5 E (q)S(q)i^t a D^(q ~ k ) ~n (S(k+)E (k) l5 S(k^)) 

Jq,k CLp 



p 2 =0 



where k± := k±p/2. Inserting the BS-equation for the vertex function in the chiral 
limit, one gets rid of the double integration in (p9|), i.e. 

2 d 



(Po^ , (So,0)P )o 



f 2 dp 2 
-1 . 



tr / 7 5 £o(fc)So(/H-)7 5 Xo(fc)£'o(£;_ 



p 2 =0 



(30) 



The last equality was obtained with the help of equation ([76]) of appendix |C]. In- 
serting (|28| ) and ([M]) into fl2~T|), we obtain the mass formula 



rz 2, 

it rat 

J TV TV 



8N r 



yj2^2\ k 2 (Z 2 — Z 2 



(2tt 



{Z 2 k 2 + Y?){Z 2 k 2 + Tti 



(31) 



This mass formula agrees with the result found earlier in reference |23[] (see also 



3.2 The Gell-Mann-Oakes-Renner relation 



Since the momentum integration at the right hand side of (pi]) is rapidly converging, 
one is tempted to conclude that the GMOR relation, i.e. 



r2 2 

ft rnl 

J TV TV 



-2m (qq) 



(32) 



does not hold in the context of the one-gluon-exchange models, because the right 
hand side of (|32D involves a divergent quark condensate whereas the left hand side is 
finite ||23|| . The apparent contradiction can be resolved by taking into account that 
the bare current mass m tends to zero for a large UV-cutoff (see (|15D) implying that 
the product of bare mass and condensate might be finite. 

The crucial quantity is the quark condensate which is defined via the trace of 
the quark propagator, i.e. 

N r r A uv ...,.„ S(A; 2 ) 



(qq) 



4n 2 



dk 2 k 



Z 2 (k 2 ) + E 2 (k 2 



(33) 



In order to show that the condensate in (|33|) is divergent in the chiral limit m^fx) = 
0, one uses the asymptotic behavior of the self-energy (|T0|) and the fact that Z(k 2 ) 
decreases to approach 1, yielding 



(qq) = (qq)(n) ln(/i 2 /A 



v QCDj 



ln(A 



rvl ^ 2 qcd) 



+ finite terms 



(34) 
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The crucial observation is that multiplying this divergent condensate with the bare 
current mass (see (^)) results in the finite value 



m(A uv ) (qq) = m R {fx) (qq)(fJ,) . 



(35) 

To establish the GMOR relation fl3~2|) , it remains to show that this finite value 
coincides up to a factor of 2 with the right hand side of the pion mass formula (BID- 



This equivalence was first observed in |£6| using the approximation Zq « Z. Below, 
we will generalize the results of [2~5 . For this purpose, we rewrite the DS-equation 
for the self-energy in the form 



£(p 2 ) = m(A uv ) + 4 / 

Jk 



d A k S(fc 2 



k *<A uv (2tt) 4 Z 2 k 2 + £ 2 



d{{p-kf 



(36) 



where d(k 2 ) = D^(k 2 )/3. For later convenience, we also introduce the DS-equation 
for the self-energy Eo in the chiral limit, i.e. 

d 4 k S (P) 



4 



d({p-kf 



(37) 



lk*<A uv (2tt) 4 Z$k 2 + E 2 , 

In order to estimate the change of the self-energy due to the presence of the bare 
mass, we subtract fl37|) from 

S(p 2 )-S (p 2 ) =m(A uv )+A 



i.e. 



d 4 k 



k 2 <A uv (27r) 



E(/c 2 



So(^ 2 



Z 2 k 2 + T 2 Z 2 k 2 + E 2 



We then multiply both sides of this equation with 

So(p 2 ) 



d((p-kf) 
(3f 



Z V)p 2 + S (p 2 ) 



and integrate over the momentum p. Using the DS-equation (|3?|), the result can be 
cast into the form 



I 



d'p ^ 2 , S(p 2 ) - S (p 2 ) 



+ 



p2<A l , v - (2tt) 4 

d A k 



S (A: 2 ) 



z V + 

E(fc 2 ) 



= m(A uv ] 
S (A: 2 ) 



d 4 p 



S (P 2 



/fc2<A av (27r) 4 

We therefore obtain the identity 



Z 2 P + S 2 Z 2 A; 2 + E 2 , 



fc 2 <A[/ v (27r) 



E (fc 2 )E(A; 2 
= m(A uv ) 



P *<A UV (2vr) 4 Z V + E 2 ^ 2 ) 

(39) 



(40) 



Z 2 *; 2 + E 2 Z 2 A; 2 + E^ 



d 4 p E (p 2 



(2tt) 4 z V + £§(p^ 
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In the product E S, one might replace the S by the self-energy S in the chiral limit, 
since the error is of higher order in the current mass Inserting then ( [4"0"D into 
the right hand side of the pion mass formula (|3lD reproduces the GMOR relation 
(§3), which completes the proof. Recently, the GMOR relation was also shown to 
hold true in effective one-gluon exchange models models which employ a running 
bare current mass |f27|. 



3.3 Numerical results 

In the table below, we finally present numerical results for the pion decay constant 
f n , pion mass m v and the quark condensate for five different values of the four quark 
coupling strength G. We also provide the results for the constituent quark mass M c , 
which is defined by the pole of the constituent quark propagator, i.e. 

p 2 = M{p 2 ) \ p2=M 2 , (41) 

where M(p 2 ) is the mass function 



GA QCD 


M c 


[-<??) 


= lGeV)] 1 / 3 


u 




1.11 


200 


MeV 


191 


MeV 


78 


MeV 


131 


MeV 


1.67 


259 


MeV 


203 


MeV 


86 


MeV 


130 


MeV 


2.22 


312 


MeV 


214 


MeV 


93 


MeV 


130 


MeV 


2.78 


362 


MeV 


223 


MeV 


99 


MeV 


130 


MeV 


3.33 


409 


MeV 


232 


MeV 


105 


MeV 


130 


MeV 



In these calculations, Aq C d ~ 500 MeV was used, and [m u + m d ](/i = 1 GeV) = 
15 MeV for up- and down-quark masses was assumed. The values in the table should 
be compared with experimental ones, i.e. 

fjfxp) = 93 MeV and m (ex P ) = 135 MeV _ 

The constituent quark mass is assumed to be approximately mjv/3 ~ 310 MeV, 
where is the nucleon mass. It is remarkable that the pion mass nearly stays 
constant when G varies. 



4 How effective is the NJL-model ? 

In this section, we establish a direct relation between our renormalizable model and 



the familiar NJL-model |[L0|| . In particular, we investigate the limits of the standard 
NJL-model Q due to its non-renormalizability. For this purpose, we will search for 
a parameter range for which our model approximately reduces to the NJL-model. 
Then we study the accuracy of this approximation in comparison to a more realistic 
choice of parameters which reproduce the phenomenological data. 
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Figure 5: The contribution of the NJL-term (dashed) and the contribution from the 
non-local term (dash-dotted) to the self-energy (full) as function of hip 2 / Aq CD for 
various values of G in units of 16ir 2 /Aj. The results are rescaled to ensure S(0) = 1. 



The effective four quark interaction of our model consists of a NJL-type contact 
interaction at low momentum transfer k < Aj and reproduces the effective perturba- 
tive one-gluon exchange at large momenta (see (0)). Consequently, the self-energy 
determined by the Dyson-Schwinger equation (|7]) receives two contributions. The 
first term coincides (up to the momentum dependence of the self-energy and to the 
presence of Z(k 2 )) with the one obtained from the NJL gap equation, while the 
second term arises from the one-gluon-exchange. The self-consistent solution to the 
equations for S(fc 2 ) and Z{k 2 ) in (0) and (§), respectively, (see figure |l|) show that 
these functions are almost constant for k < Aj, and rapidly decrease at k ph Aj to 
approach the correct asymptotic behavior of E as dictated by perturbative QCD. 
We search for a parameter range, where the second term on the right hand side of 
the DS-equation (|8|) only yields a small contribution compared with the first term. 
In this range, we can expect that observables which are insensitive to large loop 
momenta are well described by the NJL-model. The numerical analysis shows that 
such a parameter range indeed exists (see figure H). In the limit of a large effective 
four quark interaction G, the contribution to the self-energy from the non-local term 
is small. In this case, the approximately constant self-energy can be obtained from 
the reduced DS-equation i.e. 

This equation can be cast into the gap-equation of the NJL-model, i.e. 

r d 4 k M 

M = AG NJL / -— , (43) 

J\k\<A NJL (27T) 4 k 2 + M 2 
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where 



M = §M ' Gnjl = % ' Anjl = A/ ~ 500 MeV • ( 44 ) 

We therefore expect that observables involving rapidly converging loop integrals 
can be well estimated within the NJL-model, since the numerical value should be 
insensitive to the large momentum behavior. Such an observable is e.g. the pion 
decay constant defined in equation (|77|) of appendix y. 

On the other hand, the NJL-model must fail for observables which receive con- 
siderable contributions from large momenta of the loop integrals. Here, we mention 
the anomalous decays of the pion as an example. In ref. it was shown that the 



rate of the decay of the neutral pion into two photons, i.e. 

ml (a em \ 2 

°- = 16^ 9 > (45) 

can be calculated in the chiral limit independently of the details of the quark- 
propagator, i.e. 

dk 2 T, d ^n^ = / dC—±— = J , (46) 



1 + C(A; 2 )) 3 Jo (1 + C) 3 2 

where C(k 2 ) := S 2 (k 2 ) j (k 2 Z 2 (k 2 )) . The agreement of the theoretical prediction 
with the experimental value g exp = 0.504 is very good and supports the constituent 
quark picture of hadrons. Since the formula (^) is independent of the details of 
the model, it is straightforward to apply it to the case of the NJL-model. In this 
case, the momentum integration in (f46D is cutoff at Anjl, and C(k 2 ) = M 2 /k 2 . One 
obtains^ 

a%jl dk 2 ^ C - r dc 1 - 1 1 a?) 

NJL 

The typical value of the ratio of constituent quark mass and NJL-cutoff is of order 
M/Anjl ~ 1/2 implying that in the NJL-model the g-factor is underestimated by 
20%. 

The phenomenological data (see the table in the previous section) require an 
interaction strength G ~ 2.22 x 167r 2 /Aj. For this value, the NJL-term contributes 
82% to the full self-energy. This result suggests that hadronic observables calculated 
within a NJL-model, where Atvjl is tightly related to the interaction strength G^jl 



6 This result holds for an 0(4) invariant sharp cutoff. Different regularization prescriptions 
might lead to different values g. 
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by (U), should deviate from the corresponding results obtained by a renormalizable 
model by approximately 18%. This is indeed the typical accuracy of the predictions 
of the NJL-model 0. 

Finally, we mention that the NJL-propagator which arises as a limit of our 
renormalizable model possesses the structure 

S N j L (k) = — 1 - , with Z « 2.3 (48) 
Z If + 



(despite of the generic form of the gap equation In the standard NJL-model, 

the Z-factor is independent of the momentum and set to 1 in order to obtain the 
generic residuum of the quark propagator. On the other hand in renormalizable 
models, the Z-factor depends, in principle, on the interaction strength. Fortunately, 
our numerical calculations have shown (see figure |3|) that Z is approximately con- 
stant for sufficiently large coupling strength G, which is the parameter region where 
our model can be approximately reduced to an NJL-model. Therefore our inves- 



tigations suggest a NJL-model with a momentum cutoff (see equation (fl"l"D), .i.e. 



A-njl — A/ ~ 500 MeV. This value is in fact slightly smaller than the typical value, 
which arises from the fit of pion data within the NJL-model when pseudo-scalar 
axial-vector meson mixing is disregarded 



5 Summary and conclusions 

We have proposed an effective quark theory which can be regarded as an (one loop) 
renormalizable extension of the phenomenologically successful NJL-model. The ef- 
fective quark interaction of this model is given by a NJL-contact interaction at low 
momentum transfer (k < Aj), and is mediated by the one-gluon exchange augmented 
by the logarithmic QCD corrections at high energies, i.e. 



where 



G6(Aj~\k\) + ^p-edkl-A!) 



4na(Aj) 2 12vr 



kfj, k p 



G A? lQ{kj (33-2N f )\nkyAl CD - 

The intermediate energy scale Aj at which the low energy interaction matches the 
high energy one gluon exchange, is fixed by the choice of G. We demand A/ > Aqcd 
in order to account for the screening of the Landau pole 0. The model is (at least) 
one loop renormalizable. The numerical and analytical study of the asymptotic 
behavior of the quark self-energy S(fc 2 ) has yielded the correct anomalous dimensions 
for the current mass and quark condensate, which are know from perturbative QCD 
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calculations. A chiral phase transition from the perturbative chirally symmetric 
phase to a non-trivial chirally broken phase characterized by a non- vanishing quark 
condensate was found at sufficiently strong interactions, i.e. for 

/ \ 167T 2 

G > G c « 0.26 ± 0.011 . 

V ) A 2 QCD 

The numerical data for the quark mass function and the quark condensate for an 
interaction strength G close to the critical one, G c , are well fitted by 

M(p 2 = 0) » 

-(qq)(H = lGeV) « 

The critical exponent for the self-energy is consistent with 1/2, whereas the exponent 
of the quark condensate significantly deviates form 1/2. 

Within our model, the emergence of the pion as a Goldstone boson was studied. 
The Bethe-Salpeter equation for the pion vertex was solved using a derivative ex- 
pansion. We have analytically verified the Gell-Mann-Oakes-Renner relation fl32|), 
which relates the product of pion mass and decay constant to the product of current 
mass and condensate. 

We have also established the relation of our model with the familiar NJL-model 
thereby exhibiting the limits of the latter due to its non-renormalizability. Like the 
NJL-model, our model does not contain quark confinement. We believe, however, 
that quark confinement is not of much relevance for low-energy hadron dynamics 
(if the external momenta are much smaller than the quark thresholds). On the 
other hand, renormalizability seems to be much more important, in particular, for 
processes which are sensitive to large loop momenta as e.g. the anomalous decay of 
the neutral pion. Our model is obviously superior to the NJL-model in these cases. 

A Numerics 

1. Solving the Dyson-Schwinger equation 

In order to solve the coupled set of Dyson-Schwinger equations (|7||f), the integration 
over the angle (between the vectors p and k) and over the loop momentum k must 
be done numerically. For this purpose, we have discretized the momentum at n 
points in logarithmic steps and used a cubic spline interpolation for the self-energy 
S(A; 2 ) and the quark wave-function Z(k 2 ). The coupled set of differential equations 
is then solved by iteration. We used n = 20 and n = 30. It turned out that n = 20 
is sufficient for most of the applications. For numerical reasons, the momentum 



Q.UAq<-l. ( 



0.47 



0.041 A 



qcd y G 



G 



0.27 



19 



integration is cut off at Afj V = 3.3 x 10 6 Aq CD . By using also Afj V = 6.6 x 10 7 Aq CD , 
we have checked that our results are practically independent of the numerical cutoff. 
Since the functions £ and Z depend on (p + k) 2 , the numerical program needs an 
evaluation of this functions up to twice the numerical cutoff. Since the functions are 
only available for momenta smaller than A uv , we used 

{Z, E}(k 2 ) = {Z, E}(A 2 UV ) for k 2 > A 2 UV . (49) 

This procedure is part of the regularization scheme. In the scaling limit, the physical 
results are independent of the details described in ([49]). However, it turns out that 
the definition limits the numerical accuracy for the values of H(k 2 ), Z(k 2 ) for large 
values of k 2 . In order to estimate statistical and systematic numerical errors, we set 
up a second program, where we solved the Dyson- Schwinger equations (0,§) after 
the loop momentum k — > k — p was shifted. In the second program, one need not 



to resort to a truncation scheme like that in (J49) . 



From subsection |2.3j , we know that the self-energy S(/c 2 ) decreases like 1/k 2 up to 



logarithmic corrections. Close to the numerical cutoff, the self-energy takes values 
which are 10 -7 times smaller than the self-energy at zero momentum. In order to 
extract the slope of the logarithmic decrease (which is related to the quark conden- 
sate by flTU|)), a refinement of the numerical approach is necessary. It turned out to 
be convenient to introduce the auxiliary variable 

a(k 2 ) = (l + J?-\ E(k 2 ) . (50) 

\ QCD / 

At small values of k 2 , u{k 2 ) essentially coincides with E(/c 2 ). The advantage is that 
o{k 2 ) decreases only logarithmically at asymptotic large values of k 2 . To summarize, 
the second program solves the equivalent set of differential equations, i.e. 

a(p 2 ) = m + AGf *L ^ ^ CD+T>2 (51) 

a[p) m + (2 ^ )4 N ^ ^ + p ^i; 



d'k a{k 2 ) A 2 QCD +p 2 ^a{{k-p) 2 ) 



+ 4 



\k-p\>A: (2vr) 4 N(k 2 ) A 2 QCD + k 2 (k- p) 2 
AG r d 4 k Z(k 2 ) ( p-k (p-k) 2 -p 2 k 2 * 



4 f d 4 k Z(k 2 ) f o p-k , n (p-k) 2 -p 2 k T 



3 J\k- P \<Aj (2vr) 4 N(k 2 ) \ p 2 (k-p) 



3 7|fc- P |>A I (2tt) 4 N(k 2 ) \ p 2 (k-pf 
a((k — i 
(k — pY 



x 47 "y-f), (52) 



where 



N(k 2 ) = Z 3 (jfe a )jfe a + n ^- . (53) 
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It turns out that both programs produce practically the same results, which estab- 
lishes that our results are independent from the regularization scheme. The second 
program is, however, superior for calculating the quark condensate. 

2. The quark condensate 

The quark condensate can be extracted from the asymptotic behavior of the quark 
mass function in fllCf) . For this purpose, we exploit the fact that the function M(p 2 )p 2 
is almost a linear function of ln(lnp 2 ) for sufficiently large values of p 2 (see figure 
2). A line fit provides the quark condensate and the numerical estimate of the 
anomalous dimension D m . In order to estimate the error of the condensate, we use 
equation (|34|), which analytically relates the trace of the quark propagator to the 
quark condensate. We learn that the error d m — D m in the determination of the 
anomalous dimension implies that we miss the true value for the condensate by a 
factor 



a 2 1 d m —D n 
-'if 



In " uv 



A 2 



(54) 



Since we know the correct anomalous dimension analytically (d m = 4/9 from (TTJ) 
and (|15D), equation (|54j ) allows for an estimate of the numerical error, and gives rise 
to the error bars in figure |j. 

3. The critical coupling 

In order to extract the critical coupling G c (see equation ([T7|) ), we consider the self- 
energy (see figure ||)) and the quark condensate (see figure |j), respectively, at small 
values G — G c . We perform a mean square fit to the 6-8 points which are closest 
to G c . For this range of the coupling constant, we expect that the quantities £ and 
(qq) show the critical behavior according 

{E,(qq)} = c mq) (-r-1) * (55) 



c 



We then perform a mean square fit according (p5| ) to the self-energy and the conden- 
sate data points, respectively. This procedure, in particular, results in the critical 
exponents 5^ and 5(q q ) and in two (nearly identical) values for G c . The difference in 
the obtained values for G c in each case provides an estimate of the numerical error 
in G r , 



B Asymptotic form of the self-energy 

In order to estimate the contribution of the gluon exchange term (second line of 
(0))j we shift the momentum integration, yielding 

nv 2 ) « m + 4/ ^ ™_ >> 0(\ P - k\ - A,) , (56) 
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Z{p 2 



4 

3^ 



d A k 



Z{k) 



(2tt) 4 Z 2 k 2 + £ 2 



2p ■ k + tr 



4(p - A:) 2 , 



(57) 



(p - A;) 2 

We also employ a modified Landau approximation, i.e. 



« ((p - kf) 9 ((p - A;) 2 - A 2 -) « a (max(p 2 , A; 2 )) (max(p 2 , A; 2 ) - A 2 



(58) 



Numerical results confirm that this approximation is valid in the context under 
considerations of this section. The technical benefit is that we can now explicitly 
perform the angle integration in ([56]) and (|57D . In particular, the angle-integral in 
(|57|) corresponds to that of QED in Landau gauge, which is known to be zero ||18|| . It 
is therefore consistent to set Z(p 2 ) = 1 for the large values of p 2 , in agreement with 
the numerical results (see figure [I]). Within the approximation the asymptotic 
momentum dependence of the mass function M(p 2 ) coincides with that of the self- 
i.e. 



energy T,(p 2 

M{p 2 ) « S(p 2 ) , for p 2 > S 2 (0) . 
It is straightforward to evaluate the angle integration in (j56|), 



i.e. 



E(p 2 ) « m + 



a(p 2 
Tip 2 



dk 2 k 2 



E(A; 2 



+ 



A; 2 + E 2 Tr 



dk 2 a(k 2 



A; 2 + E 2 



(59) 



(60) 



where an UV-cutoff huv was introduced, that will be taken to infinity at the end. 
We will now show that the Ansatze ( |i~3D and flltf ) are self-consistent solutions to 
(^0) at large momentum p 2 . First, the case of explicit chiral symmetry breaking is 
considered. To calculate the first integral at the right hand side of (|60D, an arbitrary 
but sufficiently large scale v was introduced so that the self-energy £Qo 2 ) can be 
approximated by (|T3D) for p 2 ^> z/ 2 , i.e. 



a(p 2 
up 2 



o 



+ 



a(p 2 
up 2 

O 



k 2 



dk 2 k 



K 



a(p 2 

Tip 2 



dk 2 k 



S(A: 2 



S(A: 2 



A; 2 + S 2 
12k 



(61) 



ln a+2 p 2 /A 2 



A; 2 + S 2 

- ) + o 



(33 - 2N f 



\n a+1 p 2 /A 2 



QCD 



QCD, 



^p 2 lnp 2 /A 2 QCD/ 



The second integral on the right hand side of (|60D is 
12 



33 - 2N 



A uv du 



u ln Q+1 u/A 2 



(62) 



12 



(33 -2N f )a 



QCD 
K 



K 



\np 2 /A 2 



QCD 



\tiA 2 uv /A 2 qcd 
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We find that the term (61) is sub- leading at large p 2 . Our Ansatz solves ( |60D for 

12 d 

d m , k = m R (fi) \\n(n 2 /A 2 QCD )] ™ , (63) 



a 



33 - 2N f 
and a bare current mass 



m R (fi) \n(fi 2 /A 2 QCD ) 

m ( A UV) = dt2 /A2 • ( 64 ) 



]r\ a A 2 /A 2 

Our goal is now to show that a can be identified with the anomalous dimension of 
the current mass, and that then the high-energy behavior of the quark mass function 
as predicted by the operator product expansion (see fllOD) is obtained. 
In the absence of explicit chiral symmetry, we assume equation fllED to be the correct 
high momentum behavior of the self-energy E(p 2 ). One first estimates 

^1 f dk 2 k 2 = ^1 r dk 2 k 2 (65) 

Tip 2 Jo k 2 + S 2 up 2 Jo k 2 + S 2 



np 2 J, 2 k 2 + £ 2 m (1 - p) p 2 In? p 2 /A 2 QCD 



s p 2 \np 2 /K 2 QCD ) 

Here (3 had to be assumed smaller than 1, which will later turn out to be the case. 
The last term on the right hand side of (BUI) is 



7T 



2 

UV 



In contrast to the case of explicit chiral symmetry breaking ( p2] ) , the contribution of 
the second term (|66D is now sub-leading to the term ( |65|) (provided that < (3 < 1). 
The contribution from the upper bound Afj V in ( |5BD decreases faster than l/A^y. 
We can therefore safely remove the cutoff, i.e. Afjy — > oo, since this contribution 
decreases faster with the regulator than counter-terms. This, in particular, implies 
that we must have m R = 0, since the contribution from the regulator can not account 
for the logarithmic decrease with Ajjy of the current mass term. In order the Ansatz 
( p|) to be a solution of the Dyson-Schwinger equation in the asymptotic region (pDP, 
one must impose 

(3 = I -dm, m R (A QCD ) = . (67) 

Note that for three quark flavors Nf = 3 one finds (3 = 0.555 . . . and < (3 < 1 as 
assumed before. Our main observation here is that the Ansatz (|T6|) is the correct 
asymptotic behavior of the self-energy in the chiral limit, accompanied by the correct 
anomalous dimensions of the quark condensate as obtained by the operator product 
expansion (see (|10|)). 
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C Pion decay constant and normalization of the 
BS-vertex 



The normalization of the pion BS-vertex in (|24| ) can be derived from the observation 
that the axial- vector vertex contains a pion pole term |29] . With the assumption 



that this pole term dominates in the chiral limit, we have |15 



p 



V 

p 



2 ■ 



(6* 



where the arrows indicate the iso-spin vector. Combing this with the Ward identity 
for the axial vector vertex in the chiral limit 



s- L (k + ) l5 + l5 s-'(k.) 



(69) 



one recovers that the quark self-energy constitutes the BS-vertex in the chiral limit. 
In addition, ( |69"D provides its proper normalization, i.e. 



Po(M) 



T 75" 



So(fc) 



(70) 



A particular pion state, e.g. the ir + which we will consider in the following, is given 
by 



Po + (k) = (Pt(k) + iP 2 Q {k)) = r + l^o(k) , t + = \(t 1+ lT2 ) 



(71) 



In order to extract an explicit formula for the pion decay constant f w , we exploit the 
fact that the pion electro-magnetic charge is properly normalized to one. For this 
purpose, we study the electro-magnetic form-factor at zero photon momentum, i.e. 

2 - jf Pi(k)S(k+)r lt (k+, k + )S(k + )P (k)S(k.) (72) 

(73) 



2p fl F 7T (p 2 ,0) = tr 



- - j PS(k)S(k + )P {k)S(k-)T^, k^)S{k„ 

where k± = k ±p/2, and T fJi (k_, kJ) is the full quark photon vertex. This vertex is 
unambiguously provided by the differential form of the Ward identity, i.e. 



TJk,k) 







S-\k) . 



With the help of this equation, (ff^) can be cast into 



(74) 



2p„ 



-2tr 



I L 4W^^o( W-) + 1 1 H(k)S(k + )P (k) 



dS(k_ 
dp^ 



(75) 
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Using the decomposition S(p) = i^>cr v (p 2 ) + as(p 2 ) for the quark-propagator @, one 
realizes that both integrals at the right hand side of ( [75] ) yield the same value. One 
therefore finds 



_d_ 

dp 2 



tr / P^k)S (k + )P (k)S (k. 



p 2 =0 



(76) 



Inserting the BS- vertex (^1) in the chiral limit into (|7B|), we finally obtain an explicit 
expression for the pion decay constant, i.e. 



J 7T 



Nr. 



+ k a s a s — k oycr y + k a v a 



T 2 



-2<Js<j' s + o~ v — 2oycr y A; — k as^s 



(77) 



4_/ _/ 

V 



k 2 dk 2 



where the prime denotes the derivative of the function with respect to its argument 
k 2 . This formula was first obtained in |30| . 
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